Hydrogen atom in a magnetic field: The quadrupole moment 
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The quadrupole moment of the hydrogen atom in a mag- 
netic field B for field strengths from to 4.414 x 10 13 G is 
calculated by two different methods. The first method is vari- 
ational, and based on a single trial function. The second 
method deals with a solution of the Schrodinger equation in 
the form of a linear combination of Landau orbitals. 

PACS numbers: 31.10.+z,31.15.Pf,32.60.+i,97.20.Rp 



I. INTRODUCTION 

Plenty of works have been devoted to study of a hydro- 
gen atom in a magnetic field (see, e.g., Refs. an d 
this problem was among the first ever studied in quantum 
mechanics. To a great extent, the reason for such interest 
is due to its importance in various branches of fundamen- 
tal physics: astrophysics, spectroscopy, solid state, and 
plasma physics. From a physical point of view, the first 
appearances of the influence of a magnetic field B on 
the atom are (i) changes in binding energies, including 
the Zeeman level splitting which removes a degeneracy; 
and (ii) the development of a nonvanishing quadrupole 
moment Q a b oc B a B\, as a consequence of the defor- 
mation of the spherical-symmetrical atomic shape. In 
contrast to the former phenomenon, the latter has not 
been thoroughly studied. Meanwhile, the appearance of a 
quadrupole moment leads to a drastic change in the inter- 
action of the atoms. A standard van der Waals attraction 
which originates in the interaction of induced dipoles is 
overtaken by quadrupole-quadrupole interaction (which 
is repulsive when atoms are situated along magnetic line 
- see Refs. In many applications (for instance, 

for construction of an equation of state) , one needs to in- 
clude the effects of atom- atom interactions. For example, 
a study of pressure ionization of a strongly magnetized 
hydrogen plasma was performed in Ref. M with a sim- 
ple occupation probability model, which was based on a 
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calculation of quantum- mechanical atomic sizes |J7| . This 
model is fully adequate at sufficiently high temperatures 
T. However, in order to extend the domain of applicabil- 
ity to lower T, where the neutral fraction is large, electri- 
cal multipole interactions of atoms should be taken into 
account. Therefore, quadrupole-quadrupole interaction 
can be significant at certain plasma parameters. 

For various quantum-mechanical states of the H atom 
in a magnetic field, there have been accurate calculations 
of binding energies |^ p[, oscillator strengths pj, and 
photoionization rates Moreover, binding energies 

]7|, p2|p3| , geometrical sizes and oscillator strengths [j7|,[l3) , 
electric quadrupole transition probabilities |]l4| , and pho- 
toionization cross sections 15 have also been successfully 
calculated for an atom moving in a strong magnetic field 
(equivalent to an atom in crossed magnetic and electric 
fields), which is an essentially three-dimensional system. 
Despite this progress, up to now the quadrupole moment 
was not studied basically with probably a single excep- 
tion [|l6| . A goal of the present Brief Report is to carry 
out such a study for the ground state using (i) a vari- 
ational method and (ii) a method based on a solution 
of the Schrodinger equation by expansion in Landau or- 
bitals with coordinate-dependent coefficients. We explore 
a range of magnetic field strengths B from to the "rel- 
ativists" field B r ee m 2 e c 3 /(he) = 4.414 x 10 13 G. 



II. ASYMPTOTIC RESULTS 



Hereafter, we will measure lengths in units of clq 
Ti 2 /(m e e 2 ) = 0.529 177 A and energies in units of Ry 



7«o 



13.6057 eV. Assuming a constant uniform 



magnetic field directed along the z-axis, we take the 
vector potential A in the symmetric (axial) gauge: 
(A x ,Ay,A z ) = {B/2) (— y, x, 0). A natural parame- 
ter of the nonrelativistic theory is 7 = B/Bq, where 
B = m 2 e e 3 /(h 3 c) = 2.3505 x 10 9 G. The field is called 
"strong" if 7 > 1. 

Since the magnetic quantum number is equal to zero 
for the ground state, the Hamiltonian has the form 



r 4 
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p = x 



V 



(1) 



Because of the axial symmetry of the problem, the com- 
ponents Q a f3 of the quadrupole tensor obey the following 
relations (e.g., Ref. Jl7| ): 



1 



Qxy — Qyz 
Qxx — Qyy 



(2) 



In the weak- field limit, the usual perturbation theory 
gives fij| 



>zz = - 7 
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In the opposite case of an ultrastrong field, In 7 3> 1, 
when (x 2 ) <C (z 2 ), the longitudinal motion can be sep- 
arated, which gives rise to the one-dimensional model 
p8| . In the ground state, (z 2 ) is mainly determined by 
the exponential tail of the one-dimensional wave function, 
(z 2 ) ~ (2E)^ 1 , where E is the binding energy. Using the 
method of Hasegawa and Howard |19f| for an evaluation 
of E, we find 



Qz 



1 



21n(ln 7 ) 
Jhn) 2 ^ (ln 7 ) 3 
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at 7 



III. VARIATIONAL METHOD 



In order to construct an adequate variational trial func- 
tion ^0, we follow a recipe formulated in Refs [^0] p2[ . 
That is, the potential Vb = (A^oV^o should reproduce 
the Coulomb singularity at the origin, and the harmonic 
oscillator behavior at large distances. Furthermore, the 
trial function should have a correct functional expansion 
in coordinates at small and large distances from the ori- 
gin, as well as correct expansion in powers of B. Since 
the ground state wave function has no nodal surfaces in 
configuration space, we may write ^0 = where 4> is 
a smooth real function of coordinates. The asymptotic 
behavior of this function was calculated in Refs. p3|^6[ : 

= 7 p 2 /4 + O(r) (p->oo), (5) 
4> = r + 7 2 (Ar 3 + Brp 2 + Cr 2 + Dp 2 ) 

+0(7 V 5 ) (r-0), (6) 

where A,B,C, and D are known parameters. These 
expansions prompt to choose the following seven- 
parametric trial function 

*o = cxp I — [a 2 r 2 + (air 3 + a 2 p 2 r + a 3 p 3 + a^pr 2 )^ 
+ (h P 4 + b 2 p 2 r 2 ) 7 2 /16] 1/2 } (7) 



(cf. Refs. [ p.6[|22| ]), where a, ai-4, and bi-2 are variational 
parameters. One can check that the effective potential 
Vq corresponding to this trial function correctly repro- 
duces the potential in Eq. ([!]) at r — > (Coulomb regime) 
and at p — > 00 (Landau regime). Furthermore, Eq. (|?]) 
gives a correct functional form of the first corrections in 



powers B 2 to the exponential phase of the ground-state 
wave function (see Ref. |23f| ) and, even more importantly, 
the functional form of the first correction to the Landau 
phase factor cx Bp 2 at large distances (for a detailed dis- 
cussion, see Ref. 0). Thus, Eq. ([?]) takes into account 
the available information on the ground-state wave func- 
tion of Hamiltonian (\u) . 



IV. EXPANSION IN LANDAU ORBITALS 

The shape of the atom is close to a sphere at B -C Bq 
and to a cylinder at B 3> Bq. In the latter case, the 
expansion of the atomic wave function over the Landau 
functions is appropriate (e.g., Refs. |7|J|). 

If there were no Coulomb attraction, then the trans- 
verse part of the wave function could be described by a 
Landau function <I> ns (p, ip) [where (p is the polar angle in 
the (xy) plane], which satisfies the equation 
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dp 
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p 2 dip 2 
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(2n+ 1)7. 

(8) 



(e.g., Ref. p4j). Here n is the Landau quantum number 
and s is the negative of the z projection of the electron 
orbital momentum (n > 0, s > — n). The Landau func- 
tions form a complete orthogonal functional basis on the 
(xy) plane. 

When the atom does not move as a whole across the 
field, s is an exact quantum number. Thus a wave func- 
tion \P can be presented as 



(9) 



The sum in Eq. (||), if truncated at some n = N, can 
be considered as a variational trial function. The one- 
dimensional functions g n are to be found numerically. 
The minimum of the energy functional (\E r |?<|\f r ) implies 
zero functional derivatives: S^lH]^) /6g n (z) = (Vn). 
Taking into account Eq. (^|), one arrives at a system of 
coupled differential equations for the set of g n (z) and E, 



g n (z) +2^^,( 2 ) S „ ( (z) = (E + 2n 7 ) 5n (z), 



where 



V {s ](z) 



(10) 



pdp / dtp$Z.(p,<p)-Q nt ,,(p,<p). (11) 



The effective potentials [Eq. (|II|)] can be reduced to a 
finite sum of one-dimensional integrals feasible for nu- 
merical calculation . 
Using the relations 



2 



X 2 \ 1 

yi ) =r+r^±-(r 2 + +r 2 _), 



(12) 



^/jr-$ ns = \Jn + s + 1 $„, s +i - V™*n-M+i, 

where r± = pe ±lip , one can calculate the expectation 
values 

/oo 
z 2 | ffn (z)| 2 dz, (13) 

/OO 
[(2n + S + l)| 5 „(z)| 2 

- 2 V(n + l)(n +7+1) | 3 ;(z)5„ + i(^)|] dz. (14) 

and finally the quadrupole moment Q zz . 

At 7 3> 1 the first term n = dominates in the sum 
in Eq. (§). Hence Eq. @ results in (a; 2 ) = (y 2 ) rs 
(s + l)/7- It is worthwhile to note that neglecting all 
terms in Eq. (^) except the one at n = is equivalent 
to the adiabatic approximation used in early works (e.g., 
Refs. citeCV77,Hasegawa). 



V. RESULTS AND DISCUSSION 

The results of our calculations of the binding energy E 
and the quadrupole moment Q zz are presented in Table Q. 
When available, we compare the binding energy with the 
most accurate up-to-date results (9). 



The variational approach of Sec. 
gle seven-parametric function [Eq. 



based on a sin- 
gives very high 



relative accuracy in the binding energy on the order of 
10~ 7 at small magnetic fields which then falls to 10~ 2 at 
the largest studied magnetic fields. Basically, this corre- 
sponds to the same absolute accuracy (on the order of 
1(T 7 ) in the total energy for the whole explored range 
of magnetic fields. Two major parameters a and b\ are 
changed as a function of magnetic field in a very smooth 
and slow manner, from a ~ 1, b\ ~ 0.9 for 10 9 G to a ~ 3, 
&i ~ 0.99 at 10 13 G, respectively. Other parameters also 
vary smoothly and slowly. 



For the second method (Sec. IV), we retain n,n' = 
0, 1, . . . , 12 in the system of equations ( |l0| ) and solve it for 
the ground state at 7 > 1 using the algorithm described 
in Re f. J7j| . Then we calculate Q zz from Eq. @ using 
Eqs. and Q). 

In Table | we see that for the binding energy the 
method of expansion in the Landau orbitals turns out 
to be more accurate at 7 > 10, whereas the variational 
method of Sec. 



Ill 



is superior at lower field strengths. 
This is confirmed by a comparison with the results of 
Ref. [[I . We emphasize that our methods give very close 
results for the quadrupole moment (the deviation does 
not exceed 10%). This agrees with the qualitative be- 
havior found in Jl6| . 



The data in Table | can be approximated by the ex- 
pression 



£7 



7/4 



0.3392 +(1 + C 3 )7 7/4 ' 
where £ = 41n(l + 0.212 7 1/4 ). 



(15) 



This approximation reproduces the exact asymptotic be- 
havior: —Q zz ~ (hi7)~ 2 at 7 — > 00 and —Q zz ~ I7 2 
at 7 — > 0. Its deviation from the results in Table | does 
not exceed a few percent in the whole range of studied 
magnetic fields. 

Figure [l] shows \Q ZZ \ as a function of 7. Numerical 
results obtained as described in Sees. Ill (shown by dots) 
and IV (solid line) arc compared with perturbation the- 
ory of order B 2 and B 4 (lines marked '1' and '2', re- 
spectively) and with the fit [Eq. (|l5|)1 (dashed line). The 
quadrupole moment grows smoothly with magnetic field 
increase, reaching the maximum at 7 « 3 and then de- 
creases. For the strongly elongated atom at 7 ^> 1, 
the van der Waals constant can be roughly estimated 
as W ~ E(z 2 ) 3 oc (hi7)~ 4 . Thus W decreases at 7 — * 00 
at the same rate as Q 2 Z - This means that the distance 
R, where the van der Waals potential ~ W/R 6 becomes 
comparable with the quadrupole-quadrupole interaction 
potential ~ Q 2 /R 5 , tends to a finite value at 7 — > 00. 
Our results may have an important impact on the mod- 
eling of relatively cool neutron star atmospheres, whose 
spectra are being measured with the x-ray telescopes on- 
board the recently launched Chandra and XMM-Newton 
space observatories (e.g., Refs. p5lp6[). 
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FIG. 1. Absolute value of the quadrupole moment Q zz as 
a function of 7 = .5/(2.35 x 10 9 G) calculated by the vari- 
ational method (Sec. |lll| ) and the Landau orbital expansion 
approach (Sec. DVI): numerical results are compared with per- 
turbation theory [Eq. Q] [curve 1 corresponds to the first 
term in Eq. curve 2 to two terms] and the analytic fit 
[Eq. (^j|)]. At 7 — > the quadrupole moment Q zz tends to 
zero. 

TABLE I. Binding energy E and absolute value of the 
quadrupole moment Q zz at different magnetic fields B cal- 
culated (a) by the variational method and (b) by expansion 
in the Landau basis. Rounded-off data from Ref. |9J are given 
for comparison. 



B 




E (Ry) 




-Q zz 


(a.u.) 




(a) 


(b) 


Ref. § 


(a) 


(b) 


0.1 Bo 


1.095 05274 




1.095 05296 


0.0235 




10 9 G 


1.346 292 






0.2185 




So 


1.662 322 


1.63 


1.662 338 


0.4155 


0.38 


10 10 G 


2.614 73 


2.61 




0.5085 


0.48 


10 Bo 


3.4948 


3.490 


3.495 6 


0.4370 


0.447 


10 11 G 


5.713 


5.717 




0.2806 


0.290 


100 Bo 


7.5642 


7.579 


7.579 6 


0.2071 


0.217 


10 12 G 


11.87 


11.924 




0.1228 


0.1308 


1000 B 


15.23 


15.325 


15.324 9 


0.0915 


0.0981 


10 13 G 


22.5 


22.77 




0.0576 


0.0620 


13 r 


32.5 


32.92 




0.0380 


0.0406 
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